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Abstract 

We have studied by Monte Carlo simulations the thermal behaviour of a small 
(A^ = 13 particles) cluster described by a Heisenberg model, including nearest- 
neighbour ferromagnetic interactions and radial surface anisotropy, in an applied 
magnetic field. We have studied three different lattice structures: hexagonal close 
packed, face centered cubic and icosahedral. We show that the zero-field thermal 
behaviour depends not only on the value of the anisotropy constant but also on 
the lattice structure. The behaviour in an applied field, additionally depends, on 
the different orientations of the field with respect to the crystal axes. According to 
these relative orientations, hysteresis cycles show different step-like characteristics. 
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1 Introduction 



The magnetic properties of clusters have been studied intensely for about two 
decades. Surface and finite-size effects appear to play a key role for clusters 
containing up to several hundred atoms. An enhanced magnetic anisotropy at 
the surface has been observed for fine metallic particles [1] and predicted by 
theoretical calculations [2,3]. One source for this anisotropy is the difference 
between the coordination number at the surface and in the bulk, inducing a 
change in the crystal fields [2,4,5,6]. It has been experimentally found that 
this enhancement increases as the size of the particle decreases [1]. 

Magnetic measurements are usually carried out by passing size selected clus- 
ters through a non uniform magnetic field as in the Stern-Gerlach experiment. 
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In these experiments the cluster size is known but the crystal structure as well 
as the initial orientation of the crystal axes of the cluster with respect to the 
field gradient direction are not. 

A great deal of effort has been devoted to the theoretical and experimental 
study of transition metal (TM) clusters [7,8,9,10] which are known to have 
a relatively low anisotropy energy in the bulk compared to their exchange 
energy [2]. For TM clusters, it has been shown that it is the localised character 
of the 3d electrons at the surface which enhances the surface anisotropy with 
respect to the volume anisotropy [3]. On the contrary very little is known 
for rare earth (RE) clusters, where the anisotropy energy is expected to be 
higher. Experimental results on RE clusters show an anomalous behaviour 
(the broadening of the deflection proflle in a Stern-Gerlach experiment at a 
finite temperature) which is absent in the case of TM clusters [11,12]. 

Several studies which consider anisotropy terms have been performed, mostly 
concerning large clusters (A^ > 1000 atoms) cut to regular shapes (spherical 
or ellipsoidal) out of simple-cubic or spinel lattices. For these systems the 
competition between the surface anisotropy and the exchange term has been 
analysed as a function of temperature [13,14,15]. The influence of surface 
effects on the hysteresis cycle at zero temperature has also been shown [16]. 
Dimitrov and Wysin [17] have studied the zero-temperature behaviour of large 
particles with ferromagnetically coupled Heisenberg classical spins with either 
uniaxial random anisotropy or radial surface anisotropy. They have studied 
the two-(three-)dimensional simple-cubic lattice with the cluster having either 
circular or rectangular (spherical or cubic) geometries. They find a step-like 
hysteresis cycle which is identified as a surface effect. This is also observed 
in their zero temperature study of spherical particles cut out of an fee lattice 
structure [18]. To have a complete picture of these systems, it is of interest 
to investigate the thermal magnetic behaviour of small clusters, including an 
anisotropy energy term. 

In this work we study a simple classical Heisenberg model for an A" = 13 
cluster with radial surface anisotropy, at finite temperature and in an applied 
external field, in order to understand how the thermal magnetic behaviour 
of the cluster is affected by the existence of a surface anisotropy term. We 
will show how the interplay between the surface anisotropy and the crystal 
structure, as well as the relative orientation of the applied field with respect 
to the crystal axes, influence the magnetic behaviour of the cluster. 

This article is organized as follows: in Section 2 we describe the studied sys- 
tems; in Section 3 we give the technical details of our calculation, and in Sec- 
tion 4 we present the results. Finally, in Section 5, we summarise and discuss 
our results. 
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2 Description of the systems studied 

Among the small amount of information concerning RE clusters, reference [19] 
shows, by ab-initio electronic structure calculations that Gdi^ clusters adopt a 
hexagonal closed packed (hep) crystal structure. Numerical simulation results 
of cluster agregation starting from = 13 separated atoms placed in either 
Lennard- Jones or Gupta potentials yield very similar cohesion energies for 
hep, fee and icosahedral lattices [21]. These leads us to base our study on 
N — 13 atom clusters having one of these three crystal structures, since the 
absolute values of their cohesion energies are larger than those obtained for 
other known simple structures (such as simple-cubic or body-centered-cubic) . 

Since the magnetism in RE metals comes from the locahsed 4f orbitals, a 
Heisenberg hamiltonian is a reasonable model. Different approaches exist to 
the treatement of the spin (classical or quantum approach) and to the inter- 
action constants. In [19], based on the fact that the calculated Fermi energy 
of the band structure is in a continuum of the density of states, a classical 
RKKY-Heisenberg hamiltonian is considered to describe the Gdis cluster. On 
the other hand in [20] both quantum and classical approaches of a RKKY- 
Heisenberg hamiltonian are studied for different closed-packed N=13 clusters 
and for different geometries of the N=14 clusters. Concerning the thermal 
magnetic behaviour of the N=13 clusters, the results are qualitatively the 
same as in the classical approach developped in [19]. 

In this article we model the exchange interaction by a classical Heisenberg 
hamiltonian. In addition, a surface anisotropy term that acts only on the 
atoms at the surface to account for the reduction of the symmetry of the 
crystal at the surface is considered. 

The main sources of anisotropy are the crystal-field anisotropy and the mag- 
netic dipolar interaction. The latter, being a long range interaction, can in 
principle be neglected for small clusters. Choosing the correct term to model 
the surface anisotropy of a RE cluster is not an easy task as no experimental 
results are available. Nevertheless some hints can be taken from what is known 
for metallic clusters. In [1] it has been found that the surface anisotropy term 
increases with decreasing cluster size. 

When trying to apply these results to RE clusters one can notice that first 
principle calculations, which concern Gdis, [2,19] show a contraction of the 
surface layers that may induce a much larger crystal field anisotropy at the 
surface than in the bulk. As it has been discussed in previous articles [2,4,5], 
the lowering of the symmetry at the surface due to the missing neighbours 
produces a crystal field with a predominant axial term in the radial direction 
which can be modelled by adding to the hamiltonian a term of the form KsS^ 
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where ^ is the component along the radial vector. 

The classical RKKY-Heisenberg hamiltonian which is known to be a good 
model for bulk Gd [19] has been already studied numerically in a simplified 
version (ferromagnetic first-neighbour interactions and antiferromagnetic in- 
teractions among all other pairs) without anisotropy [23] . 

In this article we consider a complementary approach. Our aim is to study the 
infiuence of the surface anisotropy energy. For simplicity, as a first approach, 
we limit our model to the case of ferromagnetic first- neighbour couplings, 
the same at the surface and in the bulk, in order to avoid mixing the effects 
of surface anisotropy with those of competing interactions. It can be argued 
that for a very small particle, the spins on opposite faces of the surface are 
close enough to interact with each other, thus modifiying the form of the 
anisotropy. As in the case of the second neighbours interactions, this effect has 
been neglected. Since we consider clusters having an approximately spherical 
shape, we have also neglected the shape anisotropy. 

These considerations lead to the following hamiltonian which we study for the 
three crystal lattices (hep, fee, ico): 

N 

H^-J ^i-^j - Y^isi.fiif - h.Si, (1) 

<i,j> Ns i=l 

where Sj is a classical Heisenberg spin with unit length, denotes the sum 
over all nearest-neighbour pairs, Ng is the number of spins at the surface of 
the cluster (here Ng = 12), is the unit vector giving the radial direction 
from the central spin, si. J > 0, Kg > and h are the ferromagnetic interac- 
tion constant, the radial surface anisotropy constant, and the intensity of the 
applied magnetic field, measured in units oi respectively. In this article we 
fix J = 1. 



3 Calculation details 

The cluster is first simulated using standard Monte Carlo simulations for the 
temperature magnetic dependence of the system for the fee, hep and icosahe- 
dral lattices in zero field (heating and cooling simulations). 

Heating and cooling cycles of the system in constant field are simulated. The 
different considered orientations of the applied field are depicted in figure 1. 
For hep lattice these orientations are either parallel or perpendicular to the c 

— * — * 

axis of the crystal. The former is noted He — huc and the latter if„ = hun, 
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Fig. 1. Studied lattice structures, definition of Uc and Un directions of the appUed 
field for the field cooling and hysteresis simulations, (a) hep, (b)fcc, (c) icosahedral 
lattices. 



where ■u„ is a basis vector of the hexagonal central layer. As the fee lattice can 
be seen in a layered way, the same definition of the direction of the applied field 
as in the hep structure is used (see figure l.b). For the icosahedral structure, 
the Uc orientation is chosen along the direction of one of the 5-fold symmetry 
axes. The complementary direction is perpendicular to Uc and joins the 
intersection point, noted S, of the 5- fold axis with the perpendicular plane 
containing the 5 icosahedral sites, with one of these sites (see figure l.c). 
Hysteresis cycles are simulated at different temperatures, with the magnetic 
field applied in the predefined directions of each structure. 

All these simulations are done using lO^MCS/s steps for calculation after 
having discarded lO^MCS/s steps for the thcrmalisation process. Quantities 
are averaged using one configuration every lOOMCS/s in order to diminish 
statistical correlations. 

We define the average cluster magnetisation as 

N 

m^l/N <\Y,Si\>. (2) 



and the average surface magnetisation of the cluster. 



rus = l/Ns < ^ \si.ni\ >, (3) 



The energy is also calculated as a function of the temperature and the magnetic 
field. 
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Fig. 2. Ground-state energy as a function of Kg from simulated annealing for the 
hep and fee lattices. 

4 Results 

4.1 Chciraterisation of the ground states 

In order to have a better understanding of the effect of surface anisotropy at 
finite temperature and in an applied field, it is necessary to have an idea of 
the ground state of the model for the different lattice structures in zero field. 
Unfortunatelly, the analytical determination of the ground state of this hamil- 
tonian is not straightforward. The coupling between the spins and the crystal 
structure makes it impossible to perform a reduction of the number of vari- 
ables in the problem by a global rotation of the magnetic structure. Moreover, 
as this coupling is non-linear, one can not use the local field orientation search 
methods for the ground state. Wc search for the ground state numerically 
using simulated-annealing, with a decreasing power law for the temperature: 
Tj+i = (Tj)", where "i" labels the coohng step and 1.05 <a< 1.15. 

We obtain ground-state magnetisations m(0) < 1 for the three considered 
lattices due to the surface anisotropy term, which favors a canted low tem- 
perature configuration. This non-coUinear order is confirmed by the surface 
average magnetisation which saturates near ms{Q) ~ 1 for big enough values 

We note a qualitative difference between the ground-state behaviour of the hep 
and fee clusters and that of the icosahedral cluster. The ground-state energies 
and magnetisations of the hep and fee clusters are very close for low anisotropy 
and start to differ for a value around = 3, while the corresponding energy 
and magnetisation values are very different for the icosahedral cluster. Fig- 
ure 2 shows the ground-state energies for the two structures having the closest 
behaviour in the interesting region of the anisotropy, 2 < Kg < ^. 
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Figures 3 and 4 illustrate schematically the ground-state configurations found 
for the hep and fee struetures for Kg = 3, where the two lattiees start to 
differ. For the hep lattiee we find two degenerate ground states corresponding 
to two different values of the magnetic moment. For the fee cluster, only one 
magnetisation state is found. 

Two typical low temperature configurations of the hep lattice are represented 
in figures 3(a) and 3(b). In figures 3(c) and 3(d) their projections on the 
central plane may be seen (for clarity, we draw the spin orientations along 
the radial directions; the real spin configurations show slight deviations from 
these directions). The magnetic moment of the cluster (shown in figure 3(a)), 
corresponding to the ground state with the higher magnetisation, is mainly 
aligned along Uc, while for the case of the lower magnetisation ground state 
(given in figure 3(b)), the magnetic alignement is mainly along 

In figure 4(a) we show a typical ground-state configuration of the fee lattice. 
This configuration has the same magnetic moment as the hep ground state 
with the higher magnetisation (figure 3(a)). And as for this hep ground state, 
the main direction of the magnetic moment is along Uf.. 

In order to understand why the ground state is not degenerate for the fee clus- 
ter while it is degenerate for the hep cluster, we impose the low temperature 
spin configuration found for the hep lattice, shown in figure 3(b), to the fee 
lattice (respecting all the possible symmetries). This "non-physical" configu- 
ration is shown figure 4(b). Comparing figures 3(d) and 4(d), it becomes clear 
why the ground-state configuration with the lower magnetisation in the hep 
is not found in the fee cluster: the coupling between the spins in the central 
plane and those in the lower plane increases the energy in the case of the fee 
structure. 

As Kg increases, the degeneracy of the hep ground state is lifted in favour 
of the most canted ground state, i.e. with a lower magnetisation than the 
ground state of the fee cluster, and the ground states of the fee and hep 
lattices separate in energy and in magnetisation. The hep cluster is better 
able to minimize both exchange and surface energy. 

The low temperature behaviour of the icosahedral cluster always differs from 
the other two structures. The magnetic moment values m(0) of the icosahedral 
lattice at a given value of the anisotropy constant Kg are higher than those 
of the other two structures in all the anisotropy range considered. This can 
be understood by the different number of nearest neighbours at the surface. 
The sites at the surface of both hep and fee lattices have 5 nearest neighbours 
while the sites of the icosahedral lattice have 6. At fixed Kg, the ferromagnetic 
coupling favours ground-state configurations with higher magnetisations in 
the icosahedral cluster than in the other two, in agreement with our results 
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Fig. 3. Typical ground-state configurations for the hep lattiee for Kg = 3. Dotted 
lines indicate the radial directions. The real configurations show a slight deviation 
from the radial directions which we haven't drawn here for simplicity. The central 
spin is irrelevant (it follows orientation of the majority of its in-plane neighbours) 
and is not shown here for the sake of clarity, (a) Ground state corresponding to the 
case with the higher magnetic moment (mainly aligned along Uc). (b) Ground state 
corresponding to the case with the lower magnetic moment (mainly aligned along 
Un)- (c) Projection of configuration (a) on to the central plane; thick arrows indicate 
the upper plane projections and thin arrows those of the lower plane, (d) Projection 
of configuration (b) on to the central plane; upper and lower planes projections are 
coincident. 

(see figure 5). And the ground-state energies of the icosahedral lattice are 
much lower than the corresponding energies for hep and fee lattices. This is 
due to the globular geometry of the icosahedral cluster that allow a better 
compromise to minimize both exchange and surface energy. 

We will see that the similarities of both hep and fee clusters, as well as the 
specificities of the icosahedral cluster that we pointed out in this section, will 
be emphasized by the action of temperature or the application of an external 
magnetic field. 

4.-2 Zero field thermal behaviour 

We have performed heating and cooling cycles for the three considered lattices 
structures in a range of values of the anisotropy constant ^ < Kg <7 . Figure 6 
shows the average magnetisation of the particle in a cooling process in zero 
field for different Kg values and the three considered lattice structures. 
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Fig. 4. Typical ground-state configurations for the fee lattice for Kg = 3. Dotted 
lines indicate the radial directions. The real configurations show a slight deviation 
from the radial directions which we haven't drawn here for simplicity. The central 
spin is irrelevant (it follows orientation of the majority of its in-plane neighbours) 
and is not shown here for the sake of clarity (a) Ground-state configuration ; it 
has the same magnetic moment than the hep ground-state configuration (with the 
highest magnetic moment) given in figure 3(a). (b) "Non-physical configuration" ob- 
tained imposing the other hep ground-state configuration (with the lowest magnetic 
moment) given in figure 3(b) to the fee lattice ; the energy of such a configuration is 
much higher than the energy of (a), (e) Projection of configuration (a) on to the cen- 
tral plane, (d) Projection of configuration (b) on to the central plane ; thick arrows 
indicate the upper plane projections and thin arrows those of the lower plane. 

The influence of surface ansitropy on the temperature behaviour of the mag- 
netisation clearly depends on the crystal structure. For Kg = 2 the m{T) 
curves increase monotonically for the three structures. The m(T) curve of the 
icosahedral lattice is higher than those of the two other structures in all the 
measured temperature range. 

For Kg = ?) the degenerate ground state of the hep structure leads to a m(T) 
curve showing an anomalous behaviour. Which of the ground states is reached 
depends on the trajectory performed in the phase space while cooling the 
system. As the temperature is lowered, the magnetisation of the cluster grows, 
but at a value it rapidly falls down creating a peak (open circles in figure 6) . 
If the cooling process is performed very slowly, the cluster magnetisation first 
decreases and then increases again reaching the same low temperature value 
of the magnetic moment than for the fee structure (full circles in figure 6). 
Energy curves E{T) of the two runs (not shown here) superpose exactly, and 
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Fig. 5. Schematical ground-state configuration for the icosahedral lattice for Kg = 3. 
Dashed hnes indicate radial directions. For the sake of clearness, we draw the spins 
fully oriented along the radial directions in spite of the fact that the actual con- 
figurations show a slight deviation from the directions shown here. This deviation 
tends to align all the spins along the vertical C5 axis, increasing the magnetisation 
with respect to the one of the pure radial configuration. 

superpose also to the E{T) curves of the fee eluster. Obviously, for the fee 
where no degeneracy has been found, the behaviour of the m(T) curves doesn't 
show any dependence on the trajectory in the phase space during the zero- 
field cooling proccess. As expected, the two m(0) values of figure 6 for hep and 
fee lattices correspond to the magnetisation ground states found in figures 3 
and 4. 

When performing a standard MC zero-field cooling simulation of the hep lat- 
tice starting from the low magnetisation ground states found by simulated 
annealing for Kg = 3, the zero- field cooling curves are reproduced exactly and 
the peak shown in figure 6(a) appears again. 

As we have discussed in the previous section, as Kg increases, the degeneracy 
of the hep ground state is lost and the ground states of the fee and hep 
lattices separate in energy (see figure 2). For instance, for the hep cluster and 
Ks=4:, the peak is always found when cooling the system and as expected, 
the saturation value of the cluster magnetisation is smaller than for the fee 
structure (figure 6(b)). Increasing Kg diminishes the height of the peak till it 
disappears. This behaviour is a particularity of the hep structure. 

For the other lattice structures the cluster magnetisation increases monoton- 
ically when cooling and only for big values of Kg a plateau is observed. 

In all the cases the m{T) curves are ordered (from highest to lowest case) as 
follows: icosahedral, fee and hep. This fact shows that, at fixed Kg, a different 
degree of competition appears between the exchange and the surface term 
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Fig. 6. Magnetic moment as a function of temperature for hep, fee and ieosahedral 
lattice structure. Zero- field cooling simulations. Full symbols indicate a slow cooling. 

(a) Ks=2 and 3. In a slow field cooling process for the hep lattice with Ks=3, the 
magnetisation decreases when lowering T and suddenly increases again to rejoin 
the fee curve for the same Kg. The two low T states shown have the same energy. 

(b) idem, i^s=4, the slow cooling process doesn't modify the peak and the saturation 
value of m is lower than for finite T. 
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according to the considered lattice structure. We will see this fact emphasized 
by the application of an external magnetic field. 



4-3 Non-zero field thermal behaviour 

We simulated the field cooling of the system with the field applied in directions 
described in figure 1. We studied various intensities of the field ranging from 
0<h<l. 

— * 

Figure 7(a) shows m{T) for Ks='i and the two considered orientations of H. 
When the field is = hu^ {h > 0.05) the magnetisation grows monotonically 
as the temperature is decreased, and the anomaly observed in zero field disap- 
pears. On the contrary, when the field is pointing in the ■u„ direction, a peak 
is observed for the same T ^ Ta where the anomaly is found in zero field. This 
peak is confirmed by a very slow cooling process as reported in figure 7(a). 
Obviously, when the applied field is too high the peak is destroyed, and the 
magnetisation grows monotonically as T is lowered. 

The ground-state analysis of section 4.1 helps to understand why this is so. As 
decribed above, two degenerate states are found for the hep lattice for Kg = 3 
and h = 0: one corresponding to the monotonous m{T) behaviour (figure 3(a)) 
and another one corresponding to the peak in m(T) curves (figure 3(b)). The 
former has its magnetic moment mainly oriented along Uc and the latter mainly 
along Un- Then when the magnetic field is applied, one of these degenerate 
zero-field ground states is selected according to the orientation of the field. 
For the fee lattice, on the other hand, only one possibility exists (figure 4(a)) 
and the m(T) curves do not depend on the field direction. 

— * — * 

In figure 7(b) we plot \m..H\/\H\ as a function of T for different intensities 
and the two studied orientations of the applied field. It can be seen that the 
behaviour of the corresponding projection of the average magnetisation along 
the field direction depends on this direction. When the field is applied along Un 
the corresponding projection of the average magnetisation, called m„, easily 
follows the field as T is lowered. However when the field direction is Uc the 
corresponding projection, rric, grows very slowly as T is lowered saturating for 

T<Ta. 

Comparing figure 7(a) and (b) one can understand the difference in the mag- 
netic behaviour between the two chosen field directions. When the field is in 
the Uc direction, the magnetic moment grows monotonically (the peak ob- 
served in zero field disappears) while the rric projection saturates at value 
lower than one at finite T. This means that there is a non zero contribution 
to the magnetic moment which is not in the Uc direction. The situation is 
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Fig. 7. Hep lattice, Ks=3. Field-cooling simulations. He and if„ denote the magnetic 
field intensities along the directions described in figure 1(a). Full symbols curves, 
issued from slow field-cooling simulations, show that the peak is stable, (a) Magnetic 
moment as a function of temperature, (b) Projection of the magnetisation along the 
axis of the corresponding applied field. 

different when the field is in the Un direction, we can see that it is mainly m„ 
which is responsible for the peak on the ■m{T) curve. 

For the fee structure the behaviour is different. To compare with the hep 
lattice we considered the field directions as shown on the figure 1(b). The 
m(T) curves are coincident for both directions of the applied field and no 
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anomaly is observed (see figure 8(a)). This result remains true for all field 
intensities and for all values of Kg. This can be understood in terms of the 
structure of the ground state: for the fee lattice there's no degeneracy of the 
ground state so the m{T) curves do not depend on the field direction. On the 
other hand, as for the hep lattice, when the field is applied in the ■u„ direction, 
rrin follows the field more easily than does nic when the field is applied in the Uc 
direction (see figure 8(b)). The zero-temperature value of nic is lower than the 
corresponding one for m„ for all Kgi also in agreement with the hep case. This 
can also be related to the characteristics of the ground-state configurations 
in zero-field (figures 3 and 4): in both hep and fee ground states we have 
found that the spins in the central plane have a small component along the 
Uc direction, so, when a magnetic field is applied along Mc, the spins in the 
central plane find it more difficult to align with the field direction. 

The situation for the icosahedral lattice is less straightforward. First, this 
lattice has a globular structure rather than the layered structure of hep or 
fee, and so, one cannot directly identify the 5-fold symmetry with the Uc axis 
of the other two layered structures. Nevertheless, the results of field cooled 
simulations show, as for the fee case, that the orientation of the field has 
no influence in the m{T) curves. On other hand, the orientation of the fleld 
affects the projection of the magnetisation in the field direction (see figure 9). 
For this structure, it is rric which follows more easily the applied field. For 
Ks>^ nic reaches its saturation value at a finite value of T = T,. This value 
increases with Kg. This is also a sign of a non-coUinear state, confirmed by 
the high values of the surface magnetisation at zero temperature (ms(0) > 0.9 
for Ks>6). 

For a given lattice and a fixed value of the field intensity h, two different 
regimes are observed for the thermal behaviour of magnetic moment, fn{T), 
and of the projection of the magnetisation of the surface spins along the radial 
directions ms{T), according to the value of K^. For low Kg, m(T) is higher than 
ms{T) for all T. For high Kg values, the opposite behaviour is observed. There 
is then a range of K'g values where the two curves cross at a finite temperature, 
Tcross{h), leading to a non-coUinear ground state {mg{0) > m(0)). In spite of 
the fact that one cannot actually talk of a phase transition, this crossing 
indicates that the non-collinear state, where each magnetic moment has a 
strong component perpendicular to the surface, takes over from a state of large 
global magnetisation (collinear state). The range of Kg values is characteristic 
of each lattice structure. For the studied values of Kg, we have found Kg ^ 4 
for hep and fee and K^ ^ 5 for the icosahedral structure. This dependence of 
the K^ is related to the different number of nearest neighbours at the surface 
for the different structures. As already mentioned, the sites at the surface 
of both hep and fee lattices have 5 nearest neighbours while the sites of the 
icosahedral lattice have 6, which leads to a larger Kg in this latter case so as 
to counterbalance the ferromagnetic coupling. 
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Fig. 8. Fee lattice, the intensity of the applied field is h = 0.2. Field cooling sim- 
ulations. He and Hn denote the applied magnetic field parallel to the Uc and to 
the Un axis respectively (see figure 1(b)). (a)The magnetic moment as a function of 
temperature is not aff'ected by the direction of the applied field for all Kg values, 
(b) Projection of the magnetisation along the direction of the corresponding applied 
field. For all Kg values, the magnetisation follows the field more easily when it is 
applied along the Un axis. 



We have studied the Tcross{h) dependence for a given lattice at the corre- 
sponding values, when the field is applied along the Uc direction (highest 
symmetry axis for the three lattices). We observe that Tcross{h) shifts to low 
temperatures as the applied field increases. This simply reveals the fact that 
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Fig. 9. Icosahedral lattice, the intensity of the apphed field is h = 0.1. Field cooling 
simulations. He and denote the magnetic field applied parallel to the Uc and to 
the Un axis respectively (sec figure 1(c)). A plateau appears at T = (Tg increases 
with Ks). For high values of Kg { Kg > 5) large fluctuations are observed just before 
the plateau. 

when a field is applied, the thermal energy needed to flip from the non-collinear 
to the collinear state is lower. In the case of very high fields, Tc{h) 0, indi- 
cating that the ground state is already collinear. 

4-4 Hysteresis loops. 

Starting from a zero-field cooling state, we performed hysteresis loops with the 
field oriented along each of the directions described above. We observe that 
above the temperature corresponding to the anomaly of the hep lattice, no 
hysteresis is found. Then, to allow for comparison, we performed for all the 
lattices the hysteresis cycles at T < 0.1, where the hysteresis is found in hep 
structure. 

For the hep lattice no hysteresis loop is observed for Ks=l but it is already 
present for Ks=2. We observe that the field orientation along Uc gives a 
smoother loop than the orientation along Un for all Kg values. In Figure 10(a) 
it can be seen that small plateaux, and a larger coercive field, appear when 
the field is applied along the Un direction. 

For Kg=3 these characteristics persist (see Figure 10(b)). Hysteresis disap- 
pears only in the Uc direction for Kg— A and for both directions for Kg— 5. 
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Fig. 10. Hep lattice, hysteresis eycle, T = 0.01. He denotes the magnetic field applied 
parallel to the Uc axis, Hn denotes the magnetic field parallel to the Un axis (see 
figure 1(a)). Projection of the magnetisation along the axis of the corresponding 
applied field. The lines are only a guide to the eye. (a) Ks=2, (b) i^s=3. The cycle 
width is larger when the field is parallel to the Un axis. In this case the magnetisation 
remains almost unchanged till /i = 0, showing intermediate jumps before saturating 
on the opposite sens. When the field is parallel to Uc a monotonous increase of the 
magnetisation with the applied field can be observed. 

For the fee lattiee the dependenee of the eoereive field on the orientation of 
the applied field is the same as for hep, but the cycle width starts diminishing 
already at Kg—'i. 

For these two lattices, the magnetisation shows small plateaux as a function 
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of the applied field. Step-like hysteresis cycles have been observed for large 
clusters with surface anisotropy in [16,17,18] and for small antiferromagnctic 
clusters in [22]. They have been explained in terms of the simultaneous reversal 
of a group of spins. Here we show, additionally, that the location of these 
plateaux depends not only on the intensity of the field but also on its direction 
with respect to the crystal axes of the clusters. 

For the icosahedral lattice the situation is once again different from that of 
the layered structures. First, at the temperature where we performed the pre- 
vious field loops, the MC simulation fails to fiip individual spins over the high 
energy barriers. This different temperature scale is again due to the higher co- 
ordination of the surface sites, which increases the effect of the ferromagnetic 
coupling. We have then performed the field loops at T = 0.1 (the hysteresis 
disappears for Kg—S ai T — 0.15). Second, at a given T, the width of the 
hysteresis loops increases with Ks, which indicates that it is more difficult to 
reverse the magnetisation of the system as Kg increases, in contrast with the 
hep and fee cases. Again, this can be understood considering the peculiarities 
of this globular structure, and comparing its ground state with those found 
for hep and fee (see figures 3, 4 and 5). For a layered lattice one can see that 
there is a competition between the field and the anisotropy terms: when the 
cluster is (ahiiost) completely polarized, for exemple, in the +Uc direction, 
the anisotropy term in the hexagonal layer is far from being optimized (it is 
globally zero); when the field is then diminished in the reversal process, the 
anisotropy term "helps" the individual spin transitions. The situation is differ- 
ent for the icosahedral structure in which the radial character of the anisotropy 
term allows for a configuration which can follow the field without rising too 
much the anisotropy energy. When the field is reversed at low T, again in 
order to satisfy both terms, each individual spin must completely reverse its 
orientation along its radial direction, so the intermediate states involve a jump 
over the anisotropy barrier, which increases with Kg. This is exactly what we 
observe, the coercive field increases with Kg. We thus believe that the sin- 
gle spin flip algorithm used is not adapted for the study of low temperature 
hysteresis loops for the icosahedral cluster. 



5 Conclusions 

We have studied the magnetic thermal behaviour oi a, N — 13 cluster in 
a magnetic field. We have considered a classical Heisenberg model with fer- 
romagnetic interactions in a magnetic field, and a radial surface anisotropy 
term. This model has been applied to three different lattice structures, hep, 
fee and icosahedral, which have similar cohesion energies in the context of 
Lennard- Jones and Gupta potentials. 



18 



Our results show that, even in zero field, the crystal structure of the cluster 
plays an essential role on its thermal behaviour. In particular, the anomalies 
seen in the average cluster magnetisation curves as a function of the tem- 
perature for some surface anistropy constants seems characteristic of the hep 
structure. 

This anomalous magnetic behaviour has been observed on Gdi^ clusters [12] 
which are found to present a hep structure [19]. Experimental measurements 
of the magnetic moment per atom of Gdis clusters give values well below the 
one predicted by electronic structure calculations for the ground state of the 
Gdis [19] and show a peak in the m{T) curves [12,20]. 

It has been proposed that a non-coUinear arrangement of the atomic moments 
could be responsible for such a behaviour [12,19]. In [19] a simphfied version 
of a RKKY model, with no anisotropy term, is proposed as a way to ob- 
tain these non-collinear configurations. In fact, the oscillating character of the 
RKKY interactions has been replaced by a ferromagnetic nearest-neighbour 
interaction and an antiferromagnetic coupling of all other pairs in the clus- 
ter, thus renforcing the competition. In [23], using the same model on a hep 
lattice, the authors find a range of values for the competing ratio between 
ferromagnetic and antiferromagnetic interactions, 7, where the m{T) curves 
in zero field have a peak. These two studies use a classical approach. 

On the other hand, in [20] a classical and a quantum study of the model 
has been carried out. In this case they consider competing interactions: ferro- 
magnetic between first neighbours and antiferromagnetic only between second 
neighbours. They study N=13 atom closed packed clusters with: body cen- 
tered cubic, face centered cubic, hexagonal compact packing and icosahedral 
structures. They also study the N=14 hep (open shell) cluster with different 
geometries. For the closed packed (N=13) clusters, the results of the quan- 
tum approach are qualitatively the same as those of the classical one found 
in [19,23]. In all these studies the behaviour of the m(T) curves is the same: 
they show a peak where the difference 6m = m{Tmax) — Tn{Q) is less than 10% 
of the maximum value {m{Tmax)- 

In this paper, we show that the same qualitative behaviour may be induced 
by a simple hamiltonian including, in addition to the first neighbours ferro- 
magnetic interaction, a radial surface anisotropy term which accounts for a 
reduction of the symmetry of the crystal at the surface. Our results show that, 
for the hep structure, the one which Gd^^ is assumed to have [19], the compe- 
tition between the two energy terms leads to a peak in the m(T) curves for a 
certain range of Ks- For the same range of Ks the low temperature magneti- 
sation is decreased and a non-collinear configuration appears, in qualitative 
agreement with the experimental results [12] . The 5m value we found is of the 
same order than in the previous works. 
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In these two complementary approaches the constants of the model are over- 
estimated. In our case the values of leading to the peak are too big com- 
pared to the first neighbour interaction. In the RKKY approach (classical or 
quantum) the antiferromagnetic interaction intensity necessary to observe the 
peak is about 36% of the ferromagnetic one. This is around one order of mag- 
nitude higher than the ratio between the first and the second peaks of the 
RKKY interaction. In [23] this is additionnaly over-estimated by the fact that 
all the neighbours but the nearest are coupled antiferromagnetically with the 
same intensity. 

In real RE clusters one can expect both effects (competing RKKY interac- 
tions and surface anisotropy) to be present. They both contribute to the non- 
coUinear order and this may allow for the peak in the m(T) curves to be 
observed for more realistic values of Kg and 7 than those used in all these 
works. It is also useful to notice that for small clusters the experimental esti- 
mations of Curie temperature are much bigger than for the bulk. For instance, 
for G(ii3 it has been found > 500K when bulk Gd Curie temperature is 
rj-'buik _ 2Q3K [12]. This means that the cluster is able to maintain its magnetic 
order well beyond the temperature corresponding to the bulk. This pleads for 
a model which could be compatible with a canted structure and a stabilisation 
of magnetic order even for quite high temperatures. The antiferromagnetic sec- 
ond order interaction in competition with the first neighbours one goes in the 
sense of a lowering of the Tc. On the contrary, the anisotropy term contributes 
to stabilize a canted magnetic structure. This suggests that a more realistic 
model should include both therms in the hamiltonian. 

In a Stern-Gerlach experiment, the projection of the magnetisation of a sin- 
gle particle along the field gradient direction is measured. We have shown 
that this projection depends on the relative orientation between the appUed 
field and the crystal axes. Then, different relative orientations, experimentally 
unknown, will broaden Stern-Gerlach deflection profiles. Such broadened de- 
flection profiles have been observed for some RE clusters [12]. In this case, it 
is said that the clusters have a locked moment behaviour [11,12]. This means 
that the individual spins of the cluster are tightly coupled to the lattice by the 
crystal anisotropics, finding it more difficult to follow the applied field. This 
gives rise to broad deflection profiles. We show additionally, that the relative 
orientation of the applied field with respect to the crystal axes of the cluster 
also contribute to the broadening of the deflection profiles. 

We want to stress that, in this work, we deal with equilibrium properties of the 
clusters, so we are not describing the relaxation process taking place in Stern- 
Gerlach apparatus. This relaxation should be considered in the interpretation 
of the results of such experiments. This has been done by an intermediate 
approach, besides the superparamagnetic model and locked-moment model, 
in both semiclassical [24] and quantum version [25]. These models consist in 
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treating the cluster as a single magnetic moment coupled to the crystal by 
an uniaxial volume anisotropy term and to the applied field. The rotation 
degrees of freedom allow for the cluster's magnetic moment relaxation. This 
approach is complementary to ours. The anisotropy term is different and will 
not give the canted spin structure proposed for Gdis [19]. In addition, it does 
not describe the competition between exchange and anisotropy energies since 
it does not allow for individual relaxation of the spins. So no dependence on 
the lattice structure can be obtained within this approach. 

Hysteresis cycles at low temperature also show a dependence on the lattice 
structure and surface anistropy as well as on the direction of the applied field. 

Our results show that the effect of the lattice structure on the magnetic be- 
haviour cannot be neglected. Moreover, as the crystalline structure is not 
known for general N (first principles results are available only for very small 
particules), it is not excluded that, at the temperatures of the experiment, 
different structures could be present for a given N. Again, the different values 
of the low-temperature magnetisation could broaden the deflection proflle. 

Most of our results are related to the fact that, for small clusters, the struc- 
tural details become important. Authors studying very large clusters focus 
mainly on the shape of the cluster which is cut out of a lattice having simple 
cubic (or spinel) lattice structure. In these studies, it is shown how the surface 
anisotropy contributes to the existence of steps in the hysteresis loops. No 
anomaly in the m(T) curves has been reported (see [13]). Nevertheless, when 
N increases, the ratio of the number of spins at the surface to the total number 
of spins decreases and the second neigbours energy increases. It is then inter- 
esting to consider a model taking into account at the same time competing 
interactions and surface anisotropy to study how these properties evolve with 
N. This study is in progress. 
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